The Poisson-Boltzmann Theory for The Two-Plates Problem: Some Exact Results 
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The general solution to the nonlinear Poisson-Boltzmann equation for two parallel charged plates, 
either inside a symmetric elecrolyte, or inside a 2q:-q asymmetric electrolyte, is found in terms of 
Weierstrass elliptic functions. From this we derive some exact asymptotic results for the interaction 
between charged plates, as well as the exact form of the renormalized surface charge density. 
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A charged object inside an electrolyte attracts ions of 
opposite charge and repels ions of like charge. The total 
electrostatic potential, due to both the external charges 
and the electrolyte, is exponential damped as a function 
of distance from the charged object. This screening phe- 
nomenon, first studied by Debye and Huckel (DH) 
is the most essential property of ionized matter. The 
Poisson-Boltzmann (PB) equation associated with the 
screened potential is one of the most important equa- 
tions in science and engineering. Despite of this, our 
understanding of PB is very limited, largely due to its 
nonlinear nature. Analytic solutions are known only for 
a few cases, i.e. 1) one plate problem [1, 0|; and 2) one 
cylinder problem [4[, both in q:-q symmetric, and in 2q:-q 
asymmetric electrolytes. The two plate problem is the 
simplest toy problem for effective interactions between 
charged objects, and therefore has received constant at- 
tention ever since Debye-Hiickel. Its solution in symmet- 
ric electrolytes was expressed in terms of elliptic integrals 
in the classical monograph by Verwey and Overbeek Q ■ 
From this the asymptotic behaviors can be derrived in 
four different regimes, shown in Fig. [31 For a detailed 
account, see the review by Andclman [5|. The behaviors 
in cross-over regions however are not understood. More- 
over, the two-plate problem in asymmetric electrolytes 
has not been touched. In this work, we outline a number 
of new exact results for the two plate problem both in 
q:-q and in 2q:-q electrolytes. Detailed discussion will be 
published elsewhere @. 

Let us start with the Poisson-Boltzmann equations for 
the reduced, dimensionless potential ^ = fiqtp: 
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For the problem of two charged plates of infinite size, \1/ 
only depends on one coordinate z. The following first 
integrals of Eqs. ([TJ are easily obtained: 
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In the middle of two plates, d z ^(MP) — 0, therefore a 
is related to ^{MP) via 



a = cosh \fr {MP), q:-q; 
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If two plates are separated by infinite distance, ty(MP) = 
0, therefore a reduces to ao = 1 (q:-q) and ao = 1/2 
(2q:-q) respectively. The net interaction between two 
plates is proportional to 6 a = a — ao : 
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with Xsj = q 2 /4:ireT the Bjerrum length. We shall there- 
fore call 5a the reduced, dimensionless interaction. Note 
fnct vanishes as the separation becomes large. For de- 
tailed derivations of these results, see reference @. 
Now introducing a new function p(z) through 
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the first integrals Eqs. ([3]) are transformed into 
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where Eq. (|laj) applies to q:-q symmetric electrolytes, 
while Eq. (|lb[) to 2q:-q asymmetric electrolytes. j3 is 
the Boltzmann factor, q = 1.6 x 10~ 19 C the charge of 
an electron. The real space coordinates are measured in 
unit of the Debye length: 
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Eq. (0 is precisely the differential equation satisfied 
by the doubly periodic Weierstrass elliptic function 
p(z; (72, 173), with <?2 and (73 two of its invariants given 
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FIG. 1: The dimensionless potential Eqs. (|10[1 . with z a real 
variable. Solid: a = 10 in q:-q symmetric electrolyte. Dashed: 
a = 1 in 2q:-q asymmetric electrolyte, z = 0, 2ui are loga- 
rithmic singularities. Two dashed vertical lines are two maxi- 
mally separated plates with infinite surface charge density, for 
given a. Two solid red lines are two plates with finite surface 
charge density that gives the same pressure. 
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FIG. 2: The renormalized surface charge density t]r as a 
function of the bare surface charge density n for a charged 
plate as given by Eqs. (|15l) . Red solid line: q:-q. gn saturates 
at 4 in both directions. Blue dashed line: 2q:-q. r)n satu- 
rates at 6 and 6(2 — y3) in the positive and negative sides 
respectively. 
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where the prime in the summation means exclusion of 
the term with m = n = 0, and z is generally treated 
as a complex variable. p{z) is a meromorphic func- 
tion of z with an infinite number of second order poles 
z m ,n = 2towi + 2tiw2, which form a 2D rectangle lattice 
for a > cto- We shall choose u)\ to be real and W2 to be 
purely imaginary. One can easily check using Eq. (|9]) that 
p(z) is periodic with respect to two complex periods 2oj\ 
and 2l)2- p(z) = p(z + 2wi) = p(z + 2u; 2 ). There is one- 
to-one relation between the periods 2wi , 2^2 and two in- 
variants (72, 53, which can be easily computed using Wol- 
fram Mathcmatica 7. In our case, they are all uniquely 
determined by one parameter a, see Eq. ([8]). The readers 
are referred to the classic monograph by Whittaker and 
Waston Q for relevant details. 

Inverting Eqs. ©, we find the potential as 
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where the upper/lower case is for q:-q and 2q:-q elec- 
trolytes respectively Eq. (fTU)) is always positive and has 
logarithmic singularities at z = 0, 2wi, as illustrated in 
Fig. [I] It therefore describes the potential between two 
positively charged plates. Since the potential between 
plates must be a smooth function, these two singularities 
0, 2loi must be outside two plates. Hence for a given in- 
teraction a, there is a maximum separation between two 
plates, which is precisely the real period 2wi. The surface 
charge density yielding this interaction a at this maxi- 
mal separation is infinity, see the dashed vertical lines 
in Fig. [TJ This constitutes a direct proof that the inter- 
action between two plates remains finite as the surface 



charge density scales up to infinity, a property usually 
called "charge renormalization" ||. It has been repeat- 
edly observed in numerical studies of PB equation. 
, Taking the limit a — > ao, we find that two invariants 
Eq. © reduce to (1/12,-1/216) for both cases. There- 
fore Eq. (fTU)) reduces to 
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The result Eq. (I12a[) for q:-q electrolyte has been known 
since the time of Verwey and Overbeek 0]. The result 
Eq. (|12b[) for 2q:-q electrolyte was first discovered by An- 
drietti et al'vn 1976 but has remained largely unknown 
since then. Our analysis makes it clear that these solu- 
tions are related to one and the same elliptic function, 
with its real period approaching infinity. 

The general solutions of the one plate problem are ob- 
viously Eqs. (|12j) translated by an arbitrary constant z*, 
whose value is to be determined by the Neumann bound- 
ary condition. For our purpose, however, it is more con- 
venient to fix the solution to be Eqs. (|12p . and adjust the 
position of the plate zq to satisfy the boundary condition: 



av/dz\ Zo = -n. 
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where 77 = fiqa/tn = Idb/P-gc is the dimensionless 
surface charge density, with (.qc = ^/q(3a the Gouy- 
Chapman length. 

The far field asymptotics of Eqs. (fl"2"j) are given by 
±4e _z ,±6e _z respectively. On the other hand, in the 
linear PB theory, a plate with surface charge density rju 
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at the same location zq produces at z a potential 
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We identify the far- field asymptotics of the nonlinear the- 
ory with the linear theory Eq. (|14[) . and define the renor- 
malized (or effective) charge density r/u. Using Eq. (fTB"]) 
to eliminate zq in favor of 77, we find T)n(rf) as a function 
of the bare surface charge density 77: 
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which are illustrated in Fig. [2] In the weakly charged 
limit n <C 1, rja — > r/; in the strongly charged limit r/ 3> 1, 
r7# saturates at r]R (±00) = ±4 for symmetric electrolytes 
and t]r(oo) = 6,7?fl(— 00) = —6(2 — \/3) for 2q:-q asym- 
metric electrolytes. The practice of using linear theory 
with an renormalized, i.e. effective, surface charge den- 
sity in the far field, is usually called charge renormaliza- 
tion following the seminal work by Alexander et al. Here 
the most striking property is that rju saturates at a fi- 
nite value in the strongly charged limit. It shows that 
electrolytes are able to screen infinitely charged objects 
within finite distance. This should be regarded as one 
fundamental property of the Poisson-Boltzmann theory. 

If two plates are widely separated, L 3> £db, the po- 
tential between is approximately the sum of that of two 
isolated plates. The potential in the middle is 
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up to corrections smaller by powers of e~ L / lDB . Note 
that we have restored the physical unit for length here. 
Using Eqs. (Q| and Eq. ([5]), the net interaction between 
two plates is 
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where rju is determined by Eqs. (|15[) . This result applies 
to the whole region to the left of the vertical line Idb/L = 
1 in Fig. [31 which includes the intermediate regime as well 
as a major part of the Debye-Hiickel regime. 

Now consider two infinitely charged plates separated 
by a small distance L -C ioB- the potential between 
plates is high, so that the co-ion density is negligibly 
small. The corresponding term is exponentially small 
and can be dropped from the PB equation. The result- 
ing equation can be integrated easily. After some simple 
calculations, one obtains the pressure, in SI unit, as 
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FIG. 3: The parameter space for the two plate problem. 
There are four regimes where the asymptotic are well under- 
stood. Our solution in terms of Weierstrass elliptic functions 
reduces to the known asymptotic results in the corresponding 
regimes, and works in the cross-over regions shown shaded in 
the figure as well. Part of this plot was taken from the review 
by Andelman [jj. The blue line: two strongly charged plates 
with a variable distance. The red line: two weakly charged 
plates with a variable distance. 

This is called the Gouy- Chapman regime 

Consider the case of symmetric electrolytes. Assume 
L <C Idb,^gc The potential Eq. (fTU)) between plates is 
small and can be expanded into Taylor series around the 
mid-point up to quadratic order of Sz = z — uj\ : 

= *(wi) + ^"(wi)<fe a + 0(Sz 4 ). (19) 

Using the properties of Weierstrass elliptic functions:, we 
find #"(wi) = Va 2 - 1. For details, see reference [61. 
The Neumann boundary condition Eq. (1131) evaluated at 
Sz = L/2 then gives a as a function of ry, L. Substituting 
this back into Eq. ([5]) , we find: 
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This result is applicable in the top right quarter of Fig. G2 
which includes the ideal gas regime, and part of the DH 
regime. Similar analyses can also be carried out for the 
cases of asymmetric electrolytes. The results however 
are rather complicated, and not particularly illuminating. 
We therefore do not present them here. 

We shall not discuss the asymptotic behaviors in DH 
regime, as they are the same to all electrolytes and have 
been presented in various literatures. See. for example, 
the review article by Andelman Q. 

Now consider two infinitely charged plates. According 
to Fig. dU, the separation between plates is just the real 
period 2wi = 2u 1 (a) = T WH p (52(a), g 3 (a)), with g k (a) 
given by Eqs. ([S]). Here the subscript "WHP" stands 
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FIG. 4: Log-linear plots of L v.s. 8a for two infinitely charged 
plates. Black: q:-q symmetric electrolyte; Red: 77 = +00, 
2q:-q electrolyte, see Eq. (|2f I) . Green: r\ — -co, 2q:-q elec- 
trolyte, see Eq. (J22). 

for "Weierstrass Half period" . The function Twhp can 
be conveniently computed using Wolfram Mathematica 
7. The interaction between two 77 = ±00 plates in sym- 
metric electrolytes, or two 77 = +00 plates in 2q:-q elec- 
trolytes, are given by 

L = 2T W Hp(.92(a),ff 3 (a)). (21) 

The case of negatively charged plates in 2q:-q elec- 
trolytes is slightly more complicated. Here we refer the 
readers to reference 0] for details, and present the re- 
sult directly. The relation between L and a for a pair of 
negative infinitely charged plates is given by 

(a 2 a 3 1 \ 

L = 2T WH p — ,- -2Ci(a), 

V 3 27 108/ ^ U h 

—00 plates 2q:-q, (22) 

where Ci is the smaller of the two zeros of the func- 
tion p (^z + uj 2 ; %j — i^g) +a/6 in the period (0, 2lj{). 
The results Eqs. 122) are plotted in Fig. H 

For high but finite surface charge density, we have 

L = 2uji(a) - 2z , (23) 

where zq is the position of the left plate. The Neumann 
boundary condition in Eq. (|13[) can be used to solve for 
zq in terms of 77. Plugging the results back into Eq. (|23|) 
we find 
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which holds in the whole strongly charged regime where 
the Gouy-Chapman length £gc is much shorter than the 
Debye length Idb- This includes many examples in bio- 
logical physics. For a given surface charge density 77, the 



a — L curve corresponding to Eq. (|24[) can be obtained 




FIG. 5: The net pressure (a — cto) as a function of the sur- 
face charge density \n\ for given plate separation L — 2. All 
quantities are dimensionless. Red: charged plates in a sym- 
metric electrolyte. Black: positively charged plates in a 2q:-q 
asymmetric electrolyte. Green: negatively charged plates in 
a 2q:-q electrolyte. Note how multi-valence counter-ions can 
dramatically reduce the pressure between two plates. 

by a rigid shift of the saturation curve corresponding to 
Eq. (|21I22|) along the L axis. 

Finally for the general case of arbitrary surface charge 
density, we can numerically solve for a with given L, 77. 
For example, we plot for given value of L = 2£db, the 
pressure a as a function of the surface charge density 77 
in Fig. [5] for all three cases. Note how the pressure satu- 
rates as the surface charge density 77 increases. One can 
also see explicitly how divalence counter-ions dramati- 
cally reduces the repulsion between the plates, compar- 
ing with the case of mono-valence counter-ions. With 
slightly more efforts, wc can also plot a as a function of 
separation L for given surface charge density 77. We shall 
however not discuss this in detail here. 

The author thanks Leo Radzihovsky, Andy Lau, Hon- 
gru Ma, Erik Luijten, Michael Brenner, and Anatoly 
Kolomeisky for interesting discussions on electrolytes. 
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